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ABSTRACT 

We present Hubble Space Telescope NICMOS imaging of W3 IRS 5, a binary high-mass protostar. 
In addition to the two protostars, NICMOS images taken in the F222M and F160W filters show three 
new 2.22 /xm sources with very red colors; these sources fall within a region 5600 AU in diameter, 
and are coincident with a ~ 100 Mq, dense molecular clump. Two additional point sources are found 
within 0.4" (800 AU) of one of the high-mass protostars; these may be stellar companions or unresolved 
emission knots from an outflow. We propose that these sources constitute a nascent Trapezium system 
in the center of the W3 IRS 5 cluster containing as many as five proto-OB stars. This would be the 
first identification of a Trapezium still deeply embedded in its natal gas. 
Subject headings: ISM:individual(W3 IRS 5) — - infrared:stars — stars: formation 
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1. INTRODUCTION 

In the pursuit of obtaining a comprehensive under- 
standing of star formation in our Galaxy, the formation 
of stars with masses greater than 10 M© remains a chal- 
lenging problem. It has become well established that 
massive stars form not in isolation, but in close proxim- 
ity to others stars spanning the entire range of the initial 
mass function. Observationally, the gregarious nature 
of massive star formation has been established on three 
hierarchical levels. On the largest spatial scales, young 
massive stars are alm ost always found in clusters, with 
sizes around 0.5 pc (|Hodappl 119941 iSollins fc Megeatbl 
12001 . The more relevant scale for the formation of 
high-mass stars is that of a few thousand AU. On these 
spatial scales, massive stars are often found to be part 
of multiple systems and Trapezia. Trapezia, defined 
as non-hierarchical (and consequently, unstable) multi- 
ple system, are co mmonly found within OB-association s 
and open clusters ifAmbar ts umiadll95 4 . Sharplcss 1954). 
Multiplicity around high-mass stars continues to even 
smaller separations; spectroscopic observations of mas- 
sive stars show companions with separations of 10 AU or 
less [|Mermiiliod fc Garciall200l|) . 

It is not clear whether OB stars form in these multi- 
ple systems, or whether the multiple systems are formed 
from subsequent dynamical inte ractions in the dense cen- 
ters of young stellar clusters llBonnell fc Daviesl 119981 
iKroupa. Aaserth fc Hurlevl l200lF ! To date, the known 
examples of OB stars in multiple systems have already 
cleared their natal gas. If OB stars form with compan- 
ions at separations ranging from a few AU to 1000 AU, 
then these companions may play important role in me- 
diating accretion onto OB stars by modifying accretion 
disks and en velopes, by adding e xtra mass to the system 
l|Ketoll2003|) . or through mergers i|Bonnell fc Batef 2002'). 

For this reason, it is essential to search for com- 
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panions of forming massive stars at the highest avail- 
able angular resolution. A rare example of a bi- 
nary high-mass protost ellar system is W 3 IRS 5, at 
a distance of 1.8 kpc llmai et al. (2000). This in- 
frared source w as first resolved into a doubl e in- 
frared source by iHowell. McCar thy fc Low! i|198lD and 
iNeueebauer. Becklin fc Matthews! i 19821) The total lu- 
minosity is 2 x 10 5 L Q ((Campbell et alJll9^|). Radio 
observatio ns with the Very Large Array bvlWilson et alJ 
(p303l) and Ivan der Tak. Dianchi. fc Tuthilll 1)2005(1 show 
that each component of the infrared source is coincident 
with radio continuum source with sizes < 240 AU. These 
hypercompact HII regions are perhaps the youngest de- 
tectable stage of a massive star. The W3 I RS 5 region 
is also the source of at least two outflow s ijlmai et alJ 
120001 iWilson et alJl200l . iMegeath et"aT] l|1996D found 
that the W3 IRS 5 cluster is in the center of an em- 
bedded cluster of 80-240 low mass stars. We report here 
HST NICMOS observations of W3 IRS 5 with an angular 
resolution of 350 AU. 

2. OBSERVATIONS & DATA REDUCTION 

The observations were taken in Cycle 7 using cam- 
era 2 of NICMOS onb oard the Hubble Space Telescope 
(|Thompson et al.l fl998V A 3 x 3 map was made using 
a spiral-dither pattern with two dithers at each map po- 
sition. The multiaccum mode was used with a total in- 
tegration time of 128 sec for each image. The images 
were take in the F110W, F160W and F222M filters. The 
data were reduced with the standard STSCI pipeline. 
The frames were then registered and mosaiced using a 
custom program written in IDL. 

Photometry was obtained through a combination of 
aperture photometry and point spread function (PSF) 
fitting. The PSF was gene rated with the TinyTim pro- 
gram (Krist & Hook 1997). Using a custom IDL pro- 
gram, we fit the PSF to five point sources in the F222M 
mosaic: NIR1, NIR2, NIR2a, NIR2b and NIR4 (see Ta- 
ble 1). For NIR2, NIR2a and NIR2b, which are separated 
by less than 0.4", we fit three PSFs simultaneously. We 
also fit PSFs to five isolated point sources in the F222M 
mosaic and measured the photometric offset between the 
PSF fitting and the aperture photometry. We used aper- 
ture photometry for the remaining point sources in the 
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F222M and F160W mosaics, including the more isolated 
NIR3 and NIR5 sources. The total count rates were mea- 
sured over a 3 pixel aperture with a 5 to 10 pixel sky 
radius. To determine magnitudes we used the equation 
m = -2.5 x log(PHOTFNU x CR/F„{Vega)) + Corr, 
where CR is the count rate in DN sec" 1 , PHOTFNU = 
2.077 x 10" 6 and 5.502 x 10" 6 Jy s DN" 1 , F v (Vega) = 
1040.7 and 610.4 Jy, and "corr", the aperture correction, 
is -0.53 and -0.71 for the F160W and F222M bands, re- 
spectively. The coefficients are the Cycle 7 values taken 
from the online NICMOS data handbook, the aperture 
corrections were measured from the data. The resulting 
zero points for the count rate were 19.40 and 21.22 for 
the F222M and F160W bands. 

3. RESULTS 

In the center of the W3 IRS 5 cluster is a small group 
of red objects, the most prominent of which is the dou- 
ble near-IR source IRS 5 (Figures 1 and 2). In Fig- 
ure 1, we use lower limits to the colors of sources un- 
detected at F160W by adopting the lOcr detection limit 
of 20.56 mag. Within a region less than 3" in diame- 
ter, there are three sources with [F160] - [F222] > 4, 
and two additional sources which have colors [F160] — 
[F222] > 2. We designate these sources NIR 1-5. 
The five sources are coincident with a dense clump of 
molecular gas, with a mass of 100 and FWHM di- 
ameter of 5" or 9000 AU iTieftrunk et alJ 119951: 119971 
iRoberts. Crutcher. fc Troland I I1997J) . iTieftrunk et al l 
(1995) measure a column density of 5.8 x 10 cm 
toward this region with an 11" beam; the extinction 
through this c lump in the F222M band is 60 magnitudes 
l)Draineill989|) . Consequently, we find it highly unlikely 
that the detected stars are background stars. 

Images of the W3 IRS 5 system in the F160W and 
F222M filter are shown in Figure 3. iWilson et all (2003) 
found that the two brightest infrared sources, NIR1 
and NIR2, are coincident with hypercompact HII re- 
gio ns D2 and B. These were detected i n the mid-IR 
bv Ivan der Tak. Dianchi. fe Tuthilll l)2005|) . NIR3 was 
also detected in the radio continuum and mid-IR by 
Ivan der Tak. Dianchi. fc Tuthilll ()2005[) . To the east, 
there are two additional sources, NIR4 and NIR5, which 
were not detected in the radio or mid-IR. In panel C, 
we show the same image with the two hypercompact 
HII r egions subtracted us ing a PSF generated with Tiny 
Tim l|Krist fc HookllT997j) . The subtracted image reveals 
two point sources near source NIR2; these are separated 
from NIR2 by 0.3" and 0.4". We designate these sources 
NIR2a and NIR2b. In between NIR1 and NIR2 we find 
an elongated nebulosity. 

In Table 1 we give the properties of the seven un- 
resolved sources in this region. The coordinates were 
shifted so that NIR1 has the same position as the D2 
source in IWilson et alJ l|2003D . Only two sources were 
marginally detected in the F160M band, NIR1 and NIR4. 
The measured magnitudes for NIR1 and NIR2 are 21.386 
and 21.392, both having 4.5 a detections. 

The detection of NIR1, NIR2 and NIR3 at radio and 
mid-IR wavelengths essentially confirms that these are 
young OB stars embedded in the molecular core. NIR4 
does not have a detected HII region, suggesting that ei- 
ther the em ission is quenched by a high accretion rate 
ljKetoll2003h . that the star is not hot enough to produce 



TABLE 1 
Near-IR Sources 



1 

2 

2a 
2b 
3 
1 

5 



R.A. 
(J2000) 
2 25 40.784 
2 25 40.680 
2 25 40.690 
2 25 40.682 
2 25 40.727 
2 25 40.896 
2 25 41.040 



Dec 
(J2000) 
62 05 52.62 
62 05 51.63 
62 05 51.91 
62 05 51.23 
62 05 49.93 
62 05 51.71 
62 05 52.10 



F222M Sep. 1 P.A. 

Mag (unc) A.U. Deg. 

12.21 (.03) 2252 37 

13.40 (.03) - 

14.78 (.07) 519 14 

16.38 (.07) 730 177 

17.94 (.10) 3176 169 

16.20 (.07) 2780 87 

17.94 (.10) 4712 80 



Associated Radio 
or Mid-IR Source 
D2 2 /Q5 3 /K7 3 /MIR1 3 
B 2 /Q3 3 /K4 3 /MIR2 3 



Q4 3 /K6 3 /MIR3 3 



Reparation (Sep.) and Position Angle (P.A.) measu red r elative 
to NIR2. Separations assume a distance of 1830 pc Hmai et alJ 
2000) 

2 Associate d 2 cm radio source using nomenclature of 
I Claussen et alJ 119941) . also detected at 1.3 and 0.7 cm by 
IWilson et alJ 12003) 

3 Associated 1.3 cm, 0.7 cm or mid-IR source using nomneclaturc 



of lvan der Tak. Dianchi. fc Tuthfll <200" 



significant Lyman continuum radiation, or that NIR4 is 
a background star - which is unlikely given the high ex- 
tinction. Using the the tentative F160W detection, NIR4 
has a color of {FimW - F2T1M) = 5.2. Assuming 
that the emission from NIR4 is from a hot photosphere 
with no contribution from circumstellar dust, implying 
an intrinsic color of (F160JU - F222M) = 0, and adopt- 
ing a reddeni ng law of A\/E(J - H) = 2.4 x A^ 175 ) 
()Drainelll989() and a distance modulus of 11.23, we de- 
rive an absolute F222M magnitude of M222 = —1-7. Us- 
ing the tabulation of if -band magnitudes vs mass fo r 
ZAMS OB stars in lHanson. Howarth. fc Contil (l997), 
this magnitude corresponds to a spectral type B0.5V 
and a mass of 11 Mq. The lack of radio emission sug- 
gests that the star has a later spectral type than B0. 5, 
and it is likely that NIR4 is a later type star undergo- 
ing pre-main sequence evolution. For example, the Her- 
big Ae/Be star BD+65 1637, which has a spectal type 
of B3 and a mass of 9.5 M , has a dereddened abso- 
lute if -band magnitude of —1.5, similar t o the of NIR4 
(|Hillenbrand. Strom. Vrba fc Keenelll992f) . The Lyman 
continuum flux of a B3 star is 3% that of a Bl star 
l)Panagial ll973). For NIR5 there is only a lower limit 
to the measured color, (F160IU - F222M) > 2.7. As- 
suming an intrinsic color of (FlfflW - F222M) = 0.025 
(i.e. equal to the H — K color with an M3 dwarf 
iBessell fc Bre tt (1988)), the resulting absolute magni- 
tudes in the F222M band is < 3.2, corresponding to 
masses > 0.13 and > 0.3 M ffl for ages of 0.1 Myr and 
1 Myr respectively ijMuench et al.ll2~002() . This source is 
probably a low mass star. 

The two remaining point sources are found within a 
projected distance of 800 AU of NIR2; given their prox- 
imity to NIR2 and their red colors, we find it likely 
that these are stellar companions to NIR2. Neither of 
these sources would be resolved in the mid-IR imaging of 
Ivan der Tak. Dianchi. fc Tuthilll l)2005l) . but they would 
be resolved in VLA radio observations of this region. We 
see no radio counterpart to the fainter source, NIR2b, 
this source has the same magnitude as NIR5 and is prob- 
ably a low mass star. NIR2a is directly between the 



NICMOS Imaging of the W3 IRS 5 Cluster 



3 



radio sources B and A. If the 2.22 /im emission from 
NIR2a is pure photospheric emission from an embedded 
star, the star would have to be quite massive; the upper 
limit of the dereddened F22 2M magnitude is equivalent 
to th at of a ZAMS 08.5 star ^Hanson. Howarth. fc Contil 
11997ft . An alternative explanation is that the two sources 
are unresolved knots o f shocked H 2 from a bipolar flow 
originating from NIR2. iWilson et alJ l)2003ft argued that 
NIR2 /radio source B was driving an outflow in a roughly 
north-south direction; as shown in Figure 3, the two fea- 
tures are relatively aligned with the flow. These two 
explanations can be tested through narrow-band images 
in the H 2 1 -> S(l) line. 

We also report the first detection of a compact neb- 
ula between the NIR1 and NIR2. The nebula may result 
from the collision of outflows from NIR1 and NIR2, or 
from the collision of an outfl ow from one s ource with 
the envelope of the other. [Wilson et all l)2003fl and 
Ivan der T ak. Dian chi. fc Tuthilll H2005TT find evidence for 
a ionized knot moving from NIR1 at 100 km s _1 which 
may be part of an outflow moving to the north-east. If 
the proposed outflow is bipolar and centered on NIR1, 
the counter-flow would be moving to th e south-west in 
the di rection of NIR2. Observations of Clausse n et all 
( 1994) showed a radio source between NIR1 and NIR2, 
Dl, which was not dete cted in subsequent VLA imag- 
ing. IWilson et al.l l|2003ft argue that these transient ra- 
dio sources are created in outflows, in this case Dl may 
tracing the outflow responsible for the observed near- 
IR nebula. We spe culate that the X-rays detected by 
IHofner et alJ EMffl toward W3 IRS5 may originate in 
this nebula. 

3.1. Discussion 

Could W3 IRS 5 be a precursor to a Trapezium still 
embedded in its natal gas? Trapezia are defined as 
non- hierarchical multip le system of three or more stars; 
lAbt fc Corballvl ((2000) proposed that as a working rule 
for identifying Trapezia, the largest projected separa- 
tion should be no more than three times the smallest 
projected separation. The separations between the five 
near-IR sources ranges between 5600 AU and 1800 AU 
and roughly satisfy this criteria. At smaller separations, 
the detection of candidate companions around NIR2 pro- 



vides the first evidence for multiplicity at separations 
< 1000 AU around high-mass protostars. Companions 
with separations ranging from 460 AU to 16 AU have 
been detected toward the O B stars compos ing the Orion 
Nebula Trapezium system jWeigielti llilSQj) . Hence, the 
spatial distribution of the primaries, and the presence of 
companions is similar to that observed in other Trapez- 
ium system. 

The maximum projected separation of the five near- 
IR sources is 5600 AU, compar ed to 10000 AU f or the 
Trapezium in the Orion Nebula. lAbt fc Corballvl {2000) 
identified fourteen likely Trapezia, the median radius to 
the furthest outlying member of these Trapezia is 40000 
AU. Hence, the projected separations in the W3 IRS 5 
system are smaller than those observed in optically vis- 
ible Trapezia. The stellar surface density of the five 
near-IR sources composing the W3 IRS 5 Trapezia is 
10000 pc~ 2 , four times the density in the surrounding 
cluster - which is 260 pc -2 for a cluster radius of 21" 
<|Megeath et al.lll996ft . 

The small size of the W3 IRS 5 system may be a func- 
tion of its youth; such systems may expand as the natal 
gas is expelled from the system by outflows and radiation. 
Maps of the molecular gas show a compact core of dense 
gas surrounding the W3 IR S 5 with a mass of ~ 100 Mc^ 
(iTieftrymk et all 1199.51 iTieftrunk Gaume fc WilsoiJ 
119981 IRoherts. Crutcher. fc Troland I 11 997ft . The 
Lyman-continuum fluxes required by the observed 
radio continuum fluxes imply spectral types of B0. 5 
for the stars a ssocia te d with NIR1, NIR2 and NIR4 
llWilson et all 120031 Ivan der Tak Dianchi. fc Tut hill 
120051 iVacca. Garmanv fc Shulll Il996j) . Since the radio 
continuum measurements suffer from a high optical 
depth, this is a l ower limit. A B0. 5 star has a mass 
around 18 M Q i|Vacca. Garmanv fc Shull 11996ft . this 
gives a lo wer limit to t h e com bined stellar mass of 
54 M Q . iMegeath eiTaD l|1996ft found that the star 
formation efficiency, i.e. the ratio to the stellar mass 
to the total stellar and gas mass is about 18% in the 
W3 IRS 5 cluster. However, in the IRS 5 system, the 
stellar mass appears to be > 30% of the total mass. 
This high star formation efficiency may allow Trapezia 
to survive the dispersal of the natal gas and emerge as 
bound systems within a larger cluster. 
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2:25:44.00 2:25:42.67 2:25:41.33 2:25:40.00 2:25:38.67 2:25:37.33 2:25:36.00 

RA (J2000) 

Fig. 1. — The distribution of red sources detected in the F222M band. The size of the asterisk depends on the color of the source, 
ranging from 1 < (F160W - F222M) < 2 for the smallest asterisks, 2 < (FW0W - F222M) < 3, 3 < (F160W - F222M) < 4 and 
(F160W — F222M) > 4 for the largest asterisks. For sources undetected in the F160W band, a magnitude of 20.56, equivalent to a 10 a 
detection, is given as an lower limit. Sources with only lower limits to their colors are marked by circles; their colors may be redder than 
indicated by the size of the asterisk. The outer box shows the extent of the field imaged with NICMOS, the inner square is the region 
displayed in Figure 3. 
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Fig. 2. — A color composite image constructed from the F110W (blue), F160W (green) and F222M (red) mosaics of the W3 IRS 5 region, 
encompassing the whole region surveyed in the NICMOS measurements. The box shows the region displayed in Figure 3. 
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Fig. 3. — F222M (2.22 fim) and F160W (1.60 fim) images of W3 IRS 5 and the neighboring red sources and nebulosities. In Panel A we 
show the F222M image using a cube root scaling. In Panel B we show the same image, with the main NIR sources marked. The asterisks 
mark the positions of the associated radio sources D2, B, A and K6. In Panel C we show the image with the NIR1 and NIR2 sources 
subtracted. An extended nebulosity between the two sources is clearly evident. Two additional point sources sources partially hidden by 
the PSF of NIR2 are marked. The ring-like pattern is a residual from the PSF subtraction. In Panel D we show the F160W image toward 
this region, with the five IR sources marked. 



